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We show that a single imperfect fluid can be used as a source to obtain the generalized McVittie 
metric as an exact solution to Einstein's equations. The mass parameter in this metric varies with 
time thanks to a mechanism based on the presence of a temperature gradient. This fully dynamical 
solution is interpreted as an accreting black hole in an expanding universe if the metric asymptotes 
to Schwarzschild-de Sitter at temporal infinity. We present a simple but instructive example for the 
mass function and briefly discuss the structure of the apparent horizons and the past singularity. 
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The McVittie metric is a controversial solution to Ein- 
stein's equations [1] modelling a black hole in a cosmolog- 
ical FLRW background, carrying now almost 80 years of 
debates about its physical interpretation. Several papers, 
e.g. [2], [.3] and [4], pointed out that this solution has the 
correct characteristics to describe the above features. On 
the other hand, some authors [5, G] have argued that the 
McVittie solution cannot describe such systems. 

The interest for solutions of this kind, able to match 
between scales of compact objects to scales concerning 
cosmological evolution, is obvious and probably for this 
very reason many attempts in generalizing the McVittie 
solution have been made during the years, as can be seen 
in [7] and references therein. Of course, one of the fea- 
tures searched for in generalizations of this solution is the 
possibility of accreting mass on the compact object [s]. 
Despite being a quite reasonable request, especially when 
one thinks of physical objects such as stars and black 
holes, the introduction of accretion has been shown to 
lead to extremely challenging difficulties [()], even in sim- 
pler cases as in Vaidya metric and perturbations thereof 
[10]. 

In this paper we present our contribution to the subject 
by showing that, abandoning completely the assumption 
of a perfect fluid, one can in fact describe accretion (or 
evaporation) of a black hole in an FLRW background. 
We do not focus on proving that a generalized McVittie 
metric actually describes such a system, since the anal- 
ysis carried out in [2] and [4] is still applicable. Rather, 
we want to show how such a solution can be constructed, 
analyze some of its properties and present what could 
be called a toy model for the black hole evolution. We 
stick to a simple model to keep computational difficulties 
at bay, while we recognize that the imperfect fluid for- 
mulation we are presenting can accommodate for a vast 
plethora of behaviors which we plan to analyze in future 
works. In short, we look at the first generalization one 
would want to find of the non accreting McVittie black 
hole, namely one that very slowly changes its mass, keep- 
ing close to the non-accreting regime at all times. 

We shall present our solutions, specifying to a matter- 
dominated cosmological model that asymptotes to de Sit- 



ter at late times, taking the mass of the black hole to 
change continuously between two constant values. 

The McVittie metric [1] is defined by the line element 
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where a = a(t) and m is a constant. It is a solution 
to Einstein's equations initially proposed to describe a 
central object in a cosmological background driven by 
a perfect fluid with homogeneous density and inhomo- 
geneous pressure. In the limit in which H = °-/a goes 
to a constant, it reduces to the Schwarzschild-de Sitter 
solution [11], and in the limit ™/2ar ^ 1 it reduces to 
a perturbed FLRW universe with zero curvature. The 
m parameter represents the contribution to the Misner- 
Sharp mass coming from the central inhoniogeneity. 

Among the main characteristics of the spacetime de- 
scribed by this metric is the presence of a spacelike inho- 
mogeneous big bang singularity at to = 2ar, which lies 
in the causal past for a > 0. The metric has a null or 
spacelike FLRW future infinity at large r and t and two 
apparent horizons which are anti-trapped surfaces. The 
one at the lower value of r is well behaved at late times in 
the presence of a positive cosmological constant and, in 
this limit, tends to the event horizon of a Schwarzschild- 
de Sitter black hole [12]. Therefore, one can say that in 
the case in which the scale factor a asymptotes de Sit- 
ter, the metric (1) describes a black hole embedded in an 
FLRW spacetime ]2]. 

To generalize the McVittie metric, we consider a time- 
varying mass for the central object, namely to = m{t) 
in (1) ]13, 14]. Several difficulties are introduced by this 
apparently small change; the first immediate one is to 
find a reasonable fluid to enter in the energy-momentum 
tensor. 

The nonzero components of the Einstein tensor for the 
generalized McVittie line element acquire extra contribu- 
tions with respect to the original case depending on to. 
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The off-diagonal term, together with the fact that G^\, 
and G^g have to be equal, puts a stringent constraint 
on the choice of fluid. If a single perfect fluid is used 
as a source for this metric, the equality of the diagonal 
terms implies that the fluid has to be comoving [15]. This 
in turn implies that the off-diagonal term is zero, and 
therefore that m has to vanish. It follows then that no 
single perfect fluid description can be used as a source 
for the generalized McVittie. 

The problem can be alleviated with the addition of a 
second perfect fluid, which is forced though to have a 
phantom equation of state (p < —p). Moreover the two 
fluids are required to be connected by a quite unnatural 
balancing equation. 

Therefore, to find a suitable single-fluid interpretation 
for the metric (1), we require more complexity and intro- 
duce heat transport and viscosity. The most general form 
for an imperfect fiuid which is compatible with thermo- 
dynamical requirements up to first order in the gradients 
is [IG] 



T^'' = {p + p) u^u'' + p g^, - T^m^H'^'W^, 
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and where it^ is the four-velocity, T = r(a;'') the fiuid 
temperature, x the heat conductivity, rj the shear viscos- 
ity and C the bulk viscosity. 

We consider a comoving fiuid that, together with the 
spherical symmetry of the metric, implies the vanishing 
of the shear viscosity. The energy-momentum tensor thus 
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The off-diagonal component of Einstein's equations, 
given by (2b) and (7b) can be solved to find the radial de- 
pendence of the temperature, obtaining a function that 
relates the evolution of the temperature profile to the 
mass, 
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where Tao{t) is an arbitrary function of time that repre- 
sents the value of the temperature at spatial infinity. It 
is interesting to note that if m = the fiuid temperature 



is just 



Too(*) 



, which is equivalent to saying that it is in 



thermal equilibrium [17]. Therefore, the time-dependent 
black hole mass may be seen as a direct consequence 
of the fiuid being out of equilibrium. It is thanks to 
this energy transfer mechanism involving heat flow that 
the fluid can still be comoving and admit nonzero off- 
diagonal components in the metric, therefore satisfying 
what Carrera and Giulini call the spatial Ricci-isotropy 
of the Einstein tensor [■]], expressed in the first part of 
equation (2c). 

Using the remaining Einstein's equations, we can ex- 
tract solutions for the fiuid energy density in terms of the 
black hole mass 
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Of course now the freedom in our choice is expressed by 
the almost arbitrary function m. The requirements over 
m for the solution (9) to be valid are (i) that m does not 
change its sign throughout the evolution, and (ii) if d 
and rh have opposite signs, m must satisfy the constraint 
2m - H{m - 2ar) > 0. 

This defines a family of solutions, each of which is fully 
determined once Too, rn and a are chosen. Note that 
these are functions of time only, and that consequently 
the radial profiles are fully determined. While Too and a 
describe characteristics of the fluid and of the metric at 
spatial infinity, m determines the behavior of the solution 
at small r, namely the evolution of the black hole and of 
the inhomogeneous part of the fluid's energy density. 

The fluid pressure is fully determined if the above func- 
tions are given, and its expression may be obtained by 
substituting the solution for m(t), found formally inte- 
grating the chosen m, back into (2c) and (7c). The pres- 
sure differs from the original McVittie case ]2] by the 
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terms from (2c) that depend on m and by the bulk vis- 
cosity term from (7c), which introduces a further devia- 
tion from the perfect fluid case. This last effect relates 
to viscous cosmology models [IS]. 

Similarly to the original McVittie case, it is convenient 
to work in a different set of coordinates in which the 
radial coordinate coincides with the "areal radius". We 
define then f by 
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which, as in the original McVittie case, defines two 
branches [2]. We choose the branch mapping f from 
m = 2ar at f = 2m to r — > oo at f — >■ oo. The other 
branch terminates on spacelike singularities both in the 
past and in the future [2] and is thus not relevant for our 
analysis. The line element may then be cast as 
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where we have introduced the simplifying notation R — 
Jl- ^ and M = 

Y r m 

In order to determine the apparent horizons we com- 
pute the extrema of the area swept by a congruence of 
light curves. Due to spherical symmetry we only need 
to focus on radial null geodesies, which satisfy ds^ = 0. 
From (11) it immediately follows that for such curves 
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Since the area of the wavefront is given by A{f,t) ~ 
AttP , the extrema of A correspond to the solutions of 
^ = 0. In principle, the full set of solutions that 
define the surfaces we are searching for is given by 
(^)+ (^)- ^ which corresponds to an eighth or- 
der equation, as opposed to the sixth order equation one 
encounters in the original McVittie case. In our notation, 
the equation reads 
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which corresponds to gu = in (11). 

Fortunately, in the branch we are using < i? < 1, 
and the problem simplifies considerably when an accret- 
ing black hole (Af > 0) in an expanding universe {H > 0) 
is considered. In this case in fact, the outgoing null rays 
corresponding to the plus sign in (13) do not admit real 
solutions and equation (13) for the apparent horizons re- 
duces to 
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A full analysis of the causal structure is under consid- 
eration [10]. In what follows, we focus on a simple toy 
model where we take the mass of the black hole to be con- 
stant both at early and at late times, and we smoothly 
interpolate between the two static mass regimes at inter- 
mediate times. We select our model choosing a form for 
the scale factor and for the mass, 

H{t) = I + , (15) 

'l t<ta; 
m{t) = <( 1 [1 + sin {ujt + to<t<ti; (16) 

t>ti, 

where w and </> are appropriate constants to match a half 
period of the sine to the two constant mass values. When 
TO ceases to be zero the energy density and tempera- 
ture acquire gradients toward the singularity where they 
themselves go to infinity. The presence of this density 
gradient in the dynamical case avoids the rather artifi- 
cial setup of the original McVittie, whose requirement of 
a homogeneous density supported by pressure gradients 
was physically difficult to justify. Conversely, the pres- 
sure, besides showing discontinuities which are a feature 
of the oversimplification introduced in this special case, 
behaves much like in the static-mass case, going to infin- 
ity at the singularity. 

In principle, the horizon equation (14) would turn out 
to be a fourth order equation once the square roots are 
eliminated. It is important to note, though, that in the 
process of squaring spurious solutions can be introduced. 
In particular, after some work on (14) one has 
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which enforces a condition for the right-hand side to be 
positive. This constraint actually eliminates one of the 
real positive roots of the fourth order equation leaving 
only two roots that we call f+ and f_ . The solutions, as 
well as the plotted trajectories of ingoing null geodesies 
obtained by numerically solving equation (12) (following 
[4]), are plotted in Figure 1. 

Note that; for an accreting black hole, the surface f_ 
is traversable. This does not change the fact that it is 
an anti-trapping surface, rather it means that new anti- 
trapping regions are appearing above it as the mass in- 
creases and the horizon moves outwards. Figure 2 shows 
this inner region in more detail. 

In the original McVittie metric with m = the space- 
like singularity at f = 2to = f* lies to the past of all 
timelike curves if a > (a Big Bang) or to the future if 
ct < (a Big Crunch) [^]. This is seen by looking at the 
sign of ^ in the limit f — > 2to. In our dynamic mass 
case, applying such a limit to equation (12), we find 
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Figure 1. Apparent horizons for the dynamical McVittie met- 
ric with H{t) and m{t) defined in equations (f5) and (f6), 
along with trajectories of radial ingoing null geodesies. f+ 
is the outer horizon and f_ is the inner horizon. The labels 
"s" and "d" refer to the static and dynamical mass horizons, 
respectively. 
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Figure 2. Lightlike trajectories of ingoing geodesies in the 
vicinity of the inner horizon. Although the geodesies do cross 
the inner apparent horizon as it moves, they never reach the 
surface f = 2m and are eventually repelled back out of the 
horizon. 

with the next leading-order terms proportional to H van- 
ishing. This is exactly the slope of the curve f» = 2m, 
meaning that light cones become tangent to this surface, 
just as in the static-mass McVittie case. The increment 
to the variation of f with respect to the position of the 
singularity can be written as 

^(f-f,) = i?f(H-il/) + i?2(Af±l) , (19) 

which for small R depends on the sign oi H — M . If it is 
positive, as is the case in our example, every null curve 
will move away from as t increases. 

As a conclusion, we have shown that a single imperfect 
fluid can be used as a source to obtain the generalized 
McVittie metric as an exact solution to Einstein's equa- 



tions, and that the mass variation can be interpreted as 
a consequence of heat flow in the radial direction within 
the fluid. We have worked out a simple example of an 
accreting black hole to reveal its main characteristics and 
its differences with respect to the static-mass case, while 
still keeping the necessary conditions for the McVittie 
metric to be interpreted as a black hole at future infinity. 
In the case of a slow accretion rate, the main character- 
istics of the McVittie metric are still present, despite the 
shifting position of the apparent horizons and of the past 
singularity. 

In the latter part of our analysis, we have actually re- 
stricted ourselves to the simplified case of slow accretion 
when compared to the rate of expansion. In fact, if one 
moves away from this limit, we can immediately see that 
new features of the spacetime may emerge. For example, 
if one crosses the limit H = 2M, some coefficients of the 
rationalized equation (14) will vanish, drastically chang- 
ing the behavior of the apparent horizons. Furthermore, 
if one crosses H = M the spacelike character of the sin- 
gularity is no longer guaranteed. We will address these 
implications in a future work [19]. 
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